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Abstract. The pursuit of the core collapse supernova explosion mechanism 
continues. While such efforts have been undertaken over the last four decades, 
it is only in the last decade that multidimensional models have been developed, 
and only in the last few years that significant components of supernova models, 
such as the neutrino transport, have been modeled with sufficient realism. We 
can now identify what are arguably the important components that must be 
included in any realistic three-dimensional supernova model. The challenge we 
now face is to include all of these in the models. 



1. Introduction 

Four hundred years after the observation of SN1604 by Galileo Galilei, the mech- 
anism for core collapse supernova explosions remains unknown. Nonetheless, 
one-, two-, and three-dimensional simulations performed thus far have shown 
that (1) neutrino transport, (2) fluid instabilities, (3) rotation, and (4) magnetic 
fields, together with proper treatments of (5) the sub- and super- nuclear den- 
sity stellar core equation of state, (6) the neutrino interactions, and (7) gravity 
will be important. This applies to the explosion mechanism and to phenomena 
associated with supernova explosions, such as neutron star kicks, supernova neu- 
trino and gravitational wave emission, and the polarization of supernova light. 
Not surprisingly, current two- and three-dimensional models have yet to include 

(1) -(4) with sufficient realism. One-dimensional spherically symmetric mod- 
els have achieved a significant level of sophistication but, by definition, cannot 
incorporate (2)-(4), except phenomenologically. Fully general relativistic spher- 
ically sym metric simulations with B oltzmann neutrino transport do not yield 
explosions (|Liebendorfer et al. Il200ll ). demonstrating that some combination of 

(2) , (3), and (4) is required to achieve this. (While the neutrino interactions 
and equation of state are important, we do not expect changes in them to lead 
to qualitative changes in the models.) It is important to note that layering the 
dimensionality and the physics will be needed to achi eve a complete understand - 
ing of the supernova mechanism and phenomenology ( Mezzacappa et al. ll200iT ) . 



The above mentioned one-dimensional simulations, for example, provided the 
first proof that (2), (3), and (4), in some combination, are necessary ingredients 
in the recipe for explosion. 

The past modeling efforts alluded to above have illuminated t hat core col 
lapse super novae may be neutrino driven, MHD driven, or both (ISvmbalist 




1984; Wilson & Mavle 1993; Her ant et al. ll994l:lB urrows. Haves, fc Frvxell Il9 9f. . 
Jank a fc Miiller ll996l : lMacFadven fc Wooslev ll999l : lFrver fc Warrenl2004l : I.Tanka I 
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2004), but uncertainties in the current models prevent us from being able to 
single out any one of these possibilities. And it may be that more than one 
possibility is realized in Nature. Nonetheless, were a supernova neutrino driven, 
magnetic fields will likely have an impact on the dynamics of the explosion. 
Similarly, were a supernova MHD driven, the neutrino transport will dictate the 
dynamics of stellar core collapse, bounce, and the postbounce evolution, which 
in turn will create the environment in which the MHD-driven explosion would 
occur. Thus, while reduction will allow us to sort out the roles of each of the 
major physics components listed above, we will not obtain a quantitative, and 
perhaps even qualitative, understanding of core collapse supernovae until all 
components and their coupling are included in the models. 



2. Neutrino Transport 

In light of the above remarks, the "neutrino heating mechanism," whereby the 
stalled supernova shock wave is reenergized by charged-current electron neu- 
trino and antineutrino absorption on the dissociation-liberated protons and 
neutrons in the postshock flow, remains integral to the supernova paradigm, 
and neutrino transport is arguably the single most importa nt component of 
any su pernova model. Moreover, it has been shown [e.g., see iJanka fc Miiller I 
( 1996)] that the neutrino heating depends sensitively on the neutrino lumi- 
nosities, spectra, and angular distributions in the "gain" (heating) region di- 
rectly below the supernova shock wave, and these in turn are determined in 
the semitransparent region around the neutrinospheres, where neutrino trans- 
port is neither diffusive nor free streaming. In this region, only a solution to 
the neutrino Boltzmann kinetic equations can capture the three quantities of 
interest with sufficient accuracy. This is especially important when it comes to 
the neutrino spectra. The neutrino heating depends quadratically on the neu- 
trino rms energies and only linearly on the neutrino luminosities and inverse 
flux factors (a measure of the isotropy of the neutrino distributions). With 
the exception of Wilson's spherically symmetric models that invoke the doubly 
diffusive neutron finger instability in the proto-neutron star to boost the neu- 
trino luminosities, no simulation to date performed wit h multifrequency neu ~ 
trino transport has yielded an explosion, a sobering fact dWilson &; Mavldll993l: 
Swestv. Lattimer. and Mvra Ill994l: iMezzacappa et al. Ill998al: iRampp fe Jankal 

2000 JBruenn. DeNisco. fc Mezzacappa 2001; liebendorfer et al. 200ll: lThompson. Burrows, fc Pinto I 
20031 : iBuras et al.ll2003h . And this list now includes both one- and two-dimensional 
simulations. 1 Moreover, without neutron fingers, whose existence is a matter 
of curr ent debate (se e the discussion below), Wilson does not obtain explosions 
either dWilson Il200l) . 



As mentioned above, detailed spherically symmetric simulations that now 
include state of the art neutrino interactions, an industry standard equation of 
state, and multiangle, multifrequency, Boltzmann neutrino transport in full gen- 
eral relativity have been performed. Efforts by several groups are now underway 
to develop Boltzmann neutrino transport for the two- and three-dimensional 



1 A recent exception to this trend was found bv lJankal i2004f) . where a weak explosion of an 11 
Mq model was reported in the context of a two-dimensional model. 



cases (ICardall fc Mezzacappa 1 1200.^1 : ICardall 1 120041 : iLivne et al. 1 I2003T ) . The 



simulations performed thus far have been confined to the Newtonian gravity, 
O(l) limit, restricted not only by Newtonian gravity but by the exclusion of the 
0(v/c) "observer corrections" (angular aberration and frequency shift) on the 
left-hand side of the Boltzmann equation. The observer corrections are critical in 
the evolution of the comoving frame neutrino distributions, and their inclusion in 
the m o dels presents a significant techni cal challenge [e.g., see lLiebendorfer et al. I 
( 2004T ) : ICardall Mezzacappa I (J2003)], made increasingly difficult by the pro- 



liferation of terms in the Boltzmann equations as we move from one- to three- 
dimensional models and from Newtonian to general relativistic models. More- 
over, the memory and CPU requirements for three-dimensional Boltzmann neu- 
trino transport are severe and will challenge any suite of algorithms and any 
computer architecture. 

And it is now an experimental fact that neutrinos have mass and, therefore, 
mix in flavor. Observations of Solar and atmospheric neutrinos, and experiments 
at LSND, indicate there may be three independent values of the difference in 
the square of the neutrino masses (5m 2 ) and four mixing angles, which can be 
accommodated by including three active and one sterile neutrino in the models. 
Recent studies have shown that neutrino mixing in both active-active and active- 
sterile channels may occur deep in the stellar core even at small values of 5m 2 
and may alter the explosion dynam ics, supernova nucleosynthesis, and terrestria l 
neutrino signatures, or all three (|Qian Fuller I Il995l : iFuller Qiani r2004\ 
Neutrino mixing is a coherent, quantum mechanical phenomenon, unlike the 
incoherent collisional phenomena included in the Boltzmann kinetic equations. 
Therefore, a more complete (quantum kinetic) treatment of neutrino transport 
in stellar cores beyond (classical) Boltzmann transport is warranted. 



3. Fluid Instabilities 

The potential role of fluid instabilities in the post-stellar-core-bounce dynam- 
ics was first articulated decades ago and explored phenomenologically for some 
time in the context of spherically symmetric models. Fortunately, two- and 
three-dimensional models that have emerged over the past decade have now al- 
lowed us to better explore t he possibilities. (1) A s has been pointed out by a 
numb e r of authors [ e.g., see iHerant et al. I ( 1994 ) ; iBurrows. Haves, h FrvxelTl 

iFrv 



(1995h: lBuras et al.l fc003^ ; iFrver W arrc 



2004)], fluid instabilities have fun- 



damentally altered the scenario under which neutrino shock reheating occurs. 
The existence of neutrino-driven convection below the supernova shock wave al- 
lows for both an explosive scenario, in which the shock wave and material behind 
it begin to move radially outward, and continued accretion, which maintains the 
neutrino luminosities sufficiently high so as to sustain the heating. In spherical 
symmetry, explosion and accretion are mutually exclusive. While earlier studies 
speculated that convection in the proto-neutron star might also be important, 
in boosting the neutrino luminosities, the most recent and complete studies in- 
dicate that convection in the proto-neutr on star may not play a major role in 
the explosion [e.g., see lBuras et alJ (J2003)]. This conclusion, of course, does not 
necessarily extend to other instabilities in the proto-neutron star, such as doubly 
diffusive instabilities, discussed below. (2) Significant progress has been made 
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on understanding the role of doubly diffusive instabilities in core collapse super- 
novae. In particular, detailed numerical experiments performed under conditions 
culled from core collapse supernova simulations suggest that neutron fingers are 
unlikely to occur in the stellar core early after bounce, during the shock re- 
heating phase, and are therefore not e xpected to aid in generating an explosion 
((Bruenn. Ralev. Sz Mezzacappa II2004I ), 2 However, these same experiments have 
led to the discovery of a new doubly diffusive instability ( "lepto-entropy fingers" ) 
that may exist during this epoch and that may act, like neutron fingers, to boost 
the neutrinosphere luminosities. Future fully two- and three-dimensional mod- 
els will be able to investigate this possibility. (3) A completely new type of 
instability has recently been discovered in the core collapse supernova context 
( Blondin. Mezzacappa. fc DeMarino 1 2003) that may aid in generating the su- 
pernova and may be the underlying mechanism responsible for the generation of 
neutron star kicks and the polarization of supernova light. The postbounce stel- 
lar core flow is best characterized as an accretion flow through a quasi-stationary 
shock. It has been shown in two- and three-dimensional hydrodynamics stud- 
ies constructed to reflect the conditions during the postbounce shock reheating 
epoch that nonspherical perturbations of the accretion shock lead to the de- 
velopment of a "stationary accretion shock instability (SASI)." The SASI is an 
/ = 1 instability with a significant I = 2 component, and results from the nonlin- 
ear feedback between vorticity introduced by the oblique supernova shock and 
nonspherical pressure waves generated by the interacting vortices trapped in the 
postshock flow that further distort th e shoc k [this feedback mechanism was dis- 
cussed in another context bv lFoglizzo I ((20021 )] . The potential ramifications of the 
SASI for the supernova mechanism and ph enomenology wer e first elaborated in 
iBlondin. Mezzacappa. fc DeMarino"! (|2003h . Recent studies ( Janka et al. Il2004h 
confirm the existence of this instability in models that include radial ray neu- 
trino transport and explore in detail its consequences for neutron star kicks. It is 
important to note that, owing to the SASI, bipolar explosions can be produced 
even in the absence of rotation. 



4. Rotation and Magnetic Fields 

Neutrino transport in core collapse supernovae had been studied extensively over 
the last four decades, culminating recently in fully general relativistic simula- 
tions with Boltzmann neutrino transport that have essentially closed the book 
on spherically symmetric models [at least in the absence of neutrino mixing; 
iLiebendorfer et al. I (|200lh ]. And a number of two-dimensional simulations have 
now been performed over the past decade to explore the dynamics of fluid in- 
stabilities in the stellar core after bounce and their impact on the supernova 
mechanism, culminating in simulations that include sophisticated radial ray 
neu trino transport th at captures a significant amount of realism in the mod- 
els (|Buras et al.ll2003T) . In contrast, very few simulations have been performed 



2 It is important to note that when the Lattimer-Swesty equation of state iLattimer fc Swestv I 
1199 If ), used in the above studies, is used by Wil son in his simulations, the criterion for neutron 

finger instability is not satisfied ( Wils on"ll2004D . Whereas, when the Wilson equation of state 

is used, it is satisfied. This discrepancy must be explored. 
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to date that include rotation, and no contemporary, sufficiently realistic simula- 
tions have been performed that include ma gnetic fields. 

Th e work of severa l grou p s [e.g., see LBuras et al. I (I2003T): lAki vama et al.l 

(|2003h : iFrver fc WarTeil (|2004l ^: iThompson. Ouataert. fc Burrows I (|2004h ] has 
shown that rotation can significantly influence stellar core collapse and the post- 
bounce dynamics in a variety of ways. (1) Centrifugal forces will slow collapse 
along the equator and, if sufficiently large, can lead to a low-density bounce. (2) 
Gravitational binding energy will be channeled differently during collapse rela- 
tive to the spherically symmetric case. (3) The neutrinospheres will be distorted, 
and the neutrino luminosities and rms energies may be noticeably changed. (4) 
The preshock accretion ram pressure along the rotation axis and the equator 
may differ significantly. And rotation (5) will alter the development of fluid 
instabilities below the shock, (6) may provide a new source of internal energy in 
the postshock flow that may augment the energy supplied by neutrino heating, 
and (7) may have a dramatic impact on the growth of magnetic fields in the 

stellar core after bounce. 

For example, the simulations by iBuras et alJ (|2003h demonstrate clearly 
that rotation can have a dramatic effect on the postbounce shock dynamics. 
In their simulation, the average shock radius was increased by approximately a 
factor of 2 between 200 and 300 ms after bounce. Moreover, owing to centrifugal 
forces in the equatorial plane and a decreased ram pressure along the rotation 
violent overturn was observed in the postshock, convective region, and 
confined to an angular region around the rotation axis. These simulations were 
performed on a 90 degree angular grid, which will not admit the I = 1 SASI 
instability mentioned above. Future simulations that include rotation and a 
180 degree grid will allow us to more thoroughly explore whether postbounce 
instabilities and rotation act in concert in the explosion mechanism. 



A gain, the work of severa l groups [e.g., seelSvmbalistv Nl984):lMiralles. Pons. <fe Urpinl 
|): lAkivama et al.l(l2003h : lDuan fc Qianl (l2004h : lThompson. Quataert. &: Burrows I 



(2004)] has shown that, much like rotation, magnetic fields can influence stellar 
core collapse and the supernova mechanism in a number of ways: (1) The de- 
velopment of significant magnetic pressure may have an impact on stellar core 
collapse and the postbounce flow. Magnetic fields may (2) alter the develop- 
ment of fluid instabilities in the stellar core after bounce, (3) provide additional 
channels for the generation of internal energy through viscous dissipation, (4) 
alter the weak interactions in the stellar core, and (5) supply a significant, and 
perhaps dominant, fraction of the luminosity powering the explosion. Arguably, 
the fundamental question is whether the magnetic fields will organize into con- 
figurations that will drive and colli mate outflows from t he ste llar c ore. 

The pioneering simulations of iLeBlanc fc Wilson I (|1970T ) and ISvmbalistv I 



(1984) were the first to explore the evolution of stellar core magnetic fields during 



core collapse and their impact on the explosion mechanism. These simulations 
exhibited the development of a magnetic bubble deep in the core owing to the 
dramatic increase in magnetic pressure close to the rotation axis as core field 
lines are dragged inward and compressed. This magnetic bubble led to buoyant, 
bipola r outflows that culminated in bipolar explosions. However, ISvmbalistv I 
(1984) concluded that inordinately large rotation and magnetic fields strengths 
were required for an explosion to develop. The magnetic fields in the stellar 
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core can be amplified in the way just described, but the growth of field strength 
through other mechanisms must also be considered. In particular, at the time 
of the early simulations by Leblanc, Wilson, and Symbalisty, the magnetorota- 
tional instability (MRI) was not yet discovered. The MRI w as first discovered 
in th e context of the differential rotation in accretion disks (|Balbus Sz HawlevI 
1991ft but may operate in a differe ntially rotating core after stellar core bounce. 
Wheeler. Meier. Wilson! (|2002f ) delineated the possibl e roles organize d mag - 
netic fields could play in the supernova mechanism, and lAkivama etaO |2003) 
were the first to propose that the MRI could be important in the core collapse 
supernova context. They argued that, with sufficient differential rotation, the 
MRI could amplify the magnetic field strengths in the stellar core exponentially 
quickly and lead to magnetohydrodynamic luminosities ~ 10 52 erg/s, rivaling 
the neutrino luminosities. 

However, it is important to note that improved progenitor models that 
include both rotation and magnetic breaking yield iron cores that rotate at only 
~ 0.1 rad/s ((Heger. Wooslev. k, Spruit ll20 04\ Without magnetic breaking, the 
core rotation rates are 1-2 orders of magnitude faster. Therefore, magnetic 
fields lead to competing outcomes: For the MRI to increase field strengths to 
magnitudes that would be dynamically significant, perhaps dominant, would 
require significant differential rotation in the core. On the other hand, magnetic 
breaking in the core tends to slow the stellar core rotation prior to collapse. 
And, of course, two- and three-dimensional progenitor models will be required 
to better determine the differential rotation of the core at the start of collapse. 

Magnetic fields will likely play an important role in supernova dynamics. 
Whether they will take center stage as the main source of luminosity powering 
the explosion remains to be seen, and generally, their precise role will require 
realistic three-dimensional simulations. None have yet been performed. 

5. Neutrino Interactions, Sub- and Super-Nuclear Density EOS 

Improvements in modeling the macrophysics of stellar core radiation magne- 
tohydrodynamics must be matched by improvements in modeling the neutrino 
weak interactions in the stellar core and the stellar core equation of state. 

Recent simulations (jHix et al. 1 12003') have demonstrated that electron cap- 
ture during stellar core collapse is dominated by capture on nuclei, not protons, 
and that a more sophisticated treatment of capture leads to dynamically signif- 
icant changes in the size of the inner homologous core (and therefore, where the 
shock is launched) and the preshock stellar core profiles through which the shock 
will move, which will affect the subsequent shock dynamics, development of core 
fluid instabilities, and nucleosynthesis. Significant progress has been made in 
the development of nuclear structure models for heavy nuclei, but a substantial 
effort will be required to extend these models to the nuclei well above mass 100 
that populate the stellar core during collapse. 

The naive picture of neutrinos interacting on free nucleons and nuclei has 
now been replaced by the more accurate picture of neutrinos interacting with 
correlated nuclei and nucleons at sub- and super-nuclear density. Calcula- 
tions of neutrino interactions on nuclei during stellar core collapse that in- 
clude correlations among nuclei, neutrino interactions on the extended struc- 
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tures present in the nuclear "pasta" phase during the transition from nuclei 
to nuclear matter, and neutrino interactions on strongly interacting, corre- 
lat ed nucleons in the proto-neutron star h ave been completed [for a review, 
see iBurrows. Reddv. and Thompson I ( 20041 1] . Nonetheless, much more remains 
to be calculated, and these calculations will be extremely challenging. For ex- 
ample, in the case of neutrino interactions in nuclear matter, pioneering efforts 
have thus far invoked the random phase approximation (RPA), which includes 
some but not all of the correlations present in the high-density proto-neutron 
star medium. 

Simulations have been performed recently that explored the sensitivity of 
supernova models to variations in the stellar core equation of state (EOS). 
Core collapse and postb ounce simulations using the Lattimer-Swesty compress- 
ible liquid dropl et EOS (lLattimer fc Swestv Ill99lf) . the Shen et al. relativistic 
mean field EO S (She n et al II199ST ). and the Hillebrandt-Wolf Hartee-Fock EOS 
( Hillebrandt W olff 1985|) yield, for example, three-flavor neutrino lu minosi- 
ties t hat differ by tens of percent between 50 and 150 ms after bounce (|Janka I 
I2004T ). which is quantitatively significant. Therefore, further development of the 
equation of state, and such comparisons, must continue. 



6. Gravity 

One should expect the gravitational fields around proto-neutron stars to deviate 
significantly from their Newtonian values. In fact, detailed comparisons of New- 
tonian and f ully general relativistic collapse ha v e bee n performed in spherical 
symmetry bv iBruenn. DeNisco. &z Mezzacappa I ( 200 ll ^ and confirm this expec- 



tation. The results from models beginning with 15 and 25 M progenitors show 
that (1) the neutrinosphere, gain, and shock radii and (2) the infall velocities 
can differ by as much as a factor of 2 between the Newtonian and general rel- 
ativistic models. In light of these results, Newtonian multidimensional models 
can only be viewed as a stepping stone to fully general relativistic, realistic 
multidimensional models. 



7. Conclusion 

Much work remains to determine the core collapse supernova explosion mecha- 
nism and to understand all of the phenomena associated with these explosions. 
Despite significant progress and the accumulation of a wealth of understanding 
over the last four decades, as well as a number of specific clues, we still cannot an- 
swer some of the most fundamental questions. And, certainly, multidimensional 
supernova modeling is only in its infancy. However, guiding observations are 
becoming increasingly sophisticated and supercomputers are becoming increas- 
ingly powerful. Consequently, we are now presented with a unique opportunity 
to respond to this challenge. 
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